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Abstract 


This  report  describes  the  incorporation  of  an  amplifier  in  a  dipole  sensor. 
The  entire  package  is  contained  in  a  semirigid  coaxial  cable  end.  This 
configuration  provides  a  flat  E-dot  field  response  from  dc  to  2  GHz.  The 
amplification  of  the  sensor  response  before  cable  losses  allows  the  correc¬ 
tion  for  the  losses  by  convolution  in  the  frequency  domain.  The  MATLAB 
algorithms  used  to  correct  measurement  system  losses  are  described. 
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1.  Introduction 


The  accurate  measurement  of  broadband  pulse  transients  requires  com¬ 
pensation  for  effects  caused  by  measurement  system  characteristics. 
Components  of  a  measurement  system  have  different  frequency  charac¬ 
teristics.  Simple  components  such  as  cables  cause  high-frequency  rolloff. 
More  complicated  components  such  as  amplifiers,  sensors,  and  antennas 
have  more  complicated  frequency  characteristics.  From  measurements 
made  on  the  measurement  system,  a  transfer  function  is  created  that 
contains  the  information  needed  to  compensate  for  the  measurement 
system  characteristics.  This  transfer  ftmction  is  deconvolved  with  a 
measured  signal  in  the  frequency  domain  to  produce  an  accurate  meas¬ 
urement  representation. 

Because  the  compensation  is  performed  in  the  frequency  domain,  a  small 
measured  transient  may  have  portions  of  its  frequency  content  decreased 
below  the  noise  floor  by  measurement  system  attenuation.  A  portion  of 
the  information  to  correct  the  measured  signal  has  been  lost,  and  when 
the  transfer  function  is  deconvolved  with  the  measured  signal,  the  noise 
replacing  the  lost  information  is  being  corrected,  along  with  the  remain¬ 
ing  content  of  the  measured  signal. 

An  inexpensive  solution  to  this  problem  is  to  amplify  the  signal  at  the 
sensor  before  most  of  the  measurement  system  attenuation  takes  place.  In 
this  report,  the  process  of  placing  an  amplifier  at  a  dipole  junction  is 
described.  The  measurements  and  calculations  to  determine  the  system 
transfer  function  and  the  application  of  the  transfer  ftmction  to  the  meas¬ 
ured  signal  are  described.  Finally,  an  electric  field  is  measured,  corrected, 
and  compared  to  the  known  field  of  a  gigahertz-transverse  electromag¬ 
netic  (G-TEM)  cell.  This  process  is  performed  for  three  identically  fabri¬ 
cated  sensors  with  lOx,  5x,  and  no  amplification  included.  The  sensors 
produced  accurate  time-domain  electric  field  measurements  with  factors 
of  eight  and  five  increased  sensitivity  compared  to  the  nonamplified 
sensor  from  dc  to  2  GHz. 


2.  Making  the  Sensor 

The  sensor  is  a  dipole  created  at  the  end  of  a  l/8-in.-diameter  solid 
copper  shield  coaxial  cable  (coax),  which  is  used  to  support  the  dipole 
and  provide  a  return  path  for  the  signal.  The  coax  center  conductor  is  one 
leg  of  the  dipole,  and  the  shield  is  the  other  leg.  An  amplifier  is  placed 
inside  the  coax  outer  conductor  at  the  end  of  the  coax.  Pin  sockets  are 
used  to  facilitate  the  changing  of  dipole  height  without  subjecting  the 
amplifier  to  fiurther  soldering. 

The  Mini  Circuits  ERA-4  and  ERA-5  amplifiers  were  used  in  this  experi¬ 
ment.  The  small  size,  incorporated  bias-T,  and  four-pin  configuration 
made  the  amplifier  ideal  for  this  experiment.  A  great  variety  of  amplifiers 
exist  that  could  be  used  for  this  purpose  with  modifications  in  specifica- 
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tions  or  configuration.  The  four  pins  on  the  amplifier  are  input,  output, 
and  two  groimds.  A  low-temperature  soldering  iron  was  used  to  solder 
the  output  pin  to  the  center  conductor  of  the  coax.  A  short  piece  of  coax 
shield  was  slotted  at  one  end  to  fit  over  the  groimd  pins  of  the  amplifier. 
Small  pin  sockets  were  soldered  to  the  other  end  of  the  slotted  shield  and 
to  the  input  pin  of  the  amplifier.  The  slotted  shield  was  then  soldered  to 
the  end  of  the  coax  and  the  slots  were  soldered  to  the  ground  pins,  with 
care  taken  not  to  overheat  the  amplifier  leads.  The  empty  coax  was  then 
filled  to  the  end  with  epoxy  to  support  the  amplifier  and  sockets.  This 
process  is  illustrated  in  figure  1.  Wires  of  the  proper  size  to  fit  in  the 
sockets  are  bent  and  trimmed  to  form  the  dipole  legs  in  the  size  desired. 

For  the  electric  field  to  be  proportional  to  the  integral  of  the  measured 
voltage  measurement,  the  dipole  dimension  should  be  less  than  1/10  of 
the  shortest  wavelength.  The  highest  frequency  expected  is  approxi¬ 
mately  0.35  divided  by  the  risetime.  If  the  amplifier  is  within  1/10  wave¬ 
length  of  the  end  of  the  dipole  at  the  highest  frequency,  there  is  no  need  to 
worry  about  impedance  mismatches,  because  the  physical  configuration 
has  negligible  voltage  change  across  its  length. 


3.  Transfer  Function 

A  transfer  fvmction  of  a  measurement  system  component  is  performed  in 
the  time  domain  by  the  application  of  a  broadband  transient  and  the 
measurement  of  the  component  effects  on  that  transient.  A  step  waveform 
produced  by  the  Tektronix  S-52  time-domain  reflectometry  (TDR)  pulse 
generator  is  used  as  the  applied  waveform,  because  the  25-ps  risetime 
provides  ample  high-frequency  content,  and  the  low-frequency  content  is 
determined  by  the  measurement  window.  To  establish  a  baseline  of  the 
applied  signal  with  minimal  perturbation,  the  pulse  is  measured  through 
a  4-m.  coax.  A  measurement  through  the  full  system  cabling  without  the 
sensor  establishes  the  system  characteristics,  lire  complete  FFT  technique 
[1]  is  used  to  provide  a  frequency-domain  representation  of  the  step 
pulse.  The  frequency-domain  representations  of  the  system  measurement 
and  the  baseline  measurement  are  now  deconvolved  with  the  use  of  an 
iterative  method  [2]. 

The  frequency  domain  relation  between  a  signal  measurement  and  the 
actual  signal  measured  is 

Y(ja)  =  X(jQ))-H{jco),  (1) 


Figure  1.  Diagram  of 
fabricated  dipole 
sensor.  (When  all 
pieces  have  been 
soldered,  the  space 
around  the  amplifier 
is  filled  with  epoxy.) 
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where  Y(jco)  is  the  actual  signal  to  be  measured  or  the  input  to  the  meas¬ 
urement  system  component,  X(ja>)  is  the  measured  signd  data  after  the 
measurement  system  component  or  the  output  of  the  component,  and 
H(jo))  is  the  transfer  function  of  the  measurement  system  component  or  a 
representation  of  the  frequency-domain  characteristics  of  the  component. 
Values  for  Y  and  X  are  foimd  in  the  time  domain  by  measurement  of  the 
applied  signal  as  accurately  as  possible  and  by  measurement  of  the  same 
applied  signal  through  the  measurement  system  component.  To  find  the 
transfer  fimction,  one  converts  the  V  and  X  data  to  the  frequency  domain 
and  uses  these  to  solve  for  H(jco)  using  the  optimal  deconvolution 
expression  [2]; 


|X0ffl)|  +A 


(2) 


where  A  is  the  optimization  parameter  to  be  varied  in  the  iteration  proc¬ 
ess,  and  *  denotes  the  complex  conjugate. 

The  preceding  process  to  find  a  transfer  function,  given  baseline  and 
component  time  domain  data  sets,  was  implemented  in  a  MATLAB 
function  CALSTEP,  shown  in  appendix  A.  To  comply  with  local  tradition, 
we  solve  for  l/H(j(o),  or  output  over  input.  The  process  is  the  same,  with 
X  and  Y  variables  exchanged.  Progress  toward  a  solution  is  determined 
with  the  use  of  a  graphic  interface  and  knowledge  of  the  frequency 
content  of  the  applied  waveform.  A  value  of  zero  for  A  will  result  in  a 
simple  division  of  X(ja))  by  Y{j(o)  to  find  l/H(jo)),  which  is  very  noisy 
when  Y(ja))  is  small  and  X(jo))  is  not,  as  in  noisy  regions.  When  A  is  zero, 
all  information  in  the  measured  data  is  present,  including  the  noise.  As  A 
is  increased,  a  filtering  effect  occurs  when  Y(jco)  is  very  small.  This  filter¬ 
ing  occurs  in  the  high  frequencies,  where  the  magnitudes  are  very  small 
for  our  pulse-type  waveform. 


The  graphic  representation  that  is  used  to  determine  progress  plots  the 
transfer  function  for  A  equal  to  zero  and  the  present  iterative  value  of  A. 
The  plot  is  windowed  vertically  to  increase  the  sensitivity  to  the  diver¬ 
gence  of  the  curves.  A  vertical  dotted  line  is  drawn  at  the  high-frequency 
limit  of  the  applied  signal.  The  deconvolution  is  considered  optimum 
when  the  curves  overlay  to  the  high-frequency  limit.  Beyond  that  point  is 
only  information  on  noise,  which  can  be  filtered.  For  signals  with  poles 
and  zeros,  the  curves  should  overlay,  except  at  the  minimums. 


A  transfer  fimction  can  be  made  for  the  system  from  the  scope  to  the 
sensor  if  the  attenuation  is  small  enough.  If  it  is  necessary  to  create  trans¬ 
fer  functions  for  individual  components,  they  can  be  combined  into  a 
total  transfer  function  through  multiplication.  The  differences  between 
single  and  combined  transfer  functions  are  negligible.  The  scope  is  as¬ 
sumed  to  have  no  frequency  characteristics.  A  sensor  usually  has  to  be 
calibrated  separately. 


The  transfer  function  for  the  measurement  system  is  applied  to  a  meas¬ 
ured  signal  with  the  use  of  the  MATLAB  function  CAL,  which  is  shown 
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in  appendix  B.  This  function  uses  the  previously  calculated  transfer 
function  and  the  measured  time-domain  signal  as  input.  The  time- 
domain  signal  is  converted  to  the  frequency  domain,  and  the  result  is 
deconvolved  with  the  same  method  used  to  create  the  transfer  function. 
The  transfer  function  being  l/H(jco)  allows  the  use  of  the  same  calcula¬ 
tion.  The  result  is  converted  to  the  time  domain  with  an  inverse  fast 
Fourier  transform  (FFT).  The  result  of  this  process  is  shown  in  figure  2. 
This  corrected  waveform  is  then  integrated  and  processed  to  produce  the 
electric  field  transient. 


4.  The  Experiment 

The  surest  way  to  determine  the  quality  of  electric  field  measurements 
provided  by  the  amplifier/dipole  combination  is  to  measure  a  known 
field.  Our  G-TEM  cell  (EMCO  model  5317)  provides  a  known  field  to 
within  ±5  dB  when  using  a  continuous  wave  (cw)  source  above  1  GHz  [3]. 
The  imcertainty  is  due  to  standing  waves  caused  by  reflection  from  the 
absorber.  In  the  time  domain,  however,  a  narrow  pulse  could  be  meas¬ 
ured  before  the  reflections  arrive.  A  side  view  of  Ihe  G-TEM  cell  is  shown 
in  figure  3.  The  G-TEM  cell  is  an  enclosed  chamber,  shaped  like  a  horn 
antenna  with  a  horizontal  center  plane  or  septum.  This  creates  a  transmis¬ 
sion  tine  of  constant  impedance  and  increasing  volume.  The  large  end  of 
the  chamber  is  enclosed,  and  the  waU  is  covered  with  electromagnetic 
absorbing  material.  The  source  voltage  is  applied  at  the  small  end  of  the 
chamber  between  the  septum  and  the  chamber  shell.  The  G-TEM  electric 
field  at  a  point  in  the  cell  could  be  determined  from  the  applied  voltage 
and  the  height  of  the  septum  if  there  are  no  discontinuities  before  that 
point.  The  field  would  follow  a  curve  to  remain  perpendicular  to  the 
conductive  surfaces  of  the  G-TEM  cell.  Rough  measurements  to  find  the 
divergence  angle  at  the  measurement  location  indicate  that  99  percent  of 
the  electric  field  is  vertical. 


Time(s)  x10-« 
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Figure  3.  Simplified 
side  view  of  G-TEM 
EMCO  model  5317 
showing  physical 
placement  of  septum 
and  cell  walls. 


G-TEM 

center 

line 


A  TDR  measurement  was  performed  to  determine  the  impedance  charac¬ 
teristics  of  the  G-TEM  cell.  The  results,  plotted  in  figure  4,  show  a  flat 
impedance  imtil  the  absorber.  The  cell  begins  at  the  input  cable  end, 
which  can  be  seen  by  the  open  circuit  reflection.  Comparison  of  the 
curves  before  the  open  circuit  indicates  that  the  spikes  are  noise,  but  the 
ripple  is  not.  A  steady  increase  in  resistance  is  indicated  by  the  increasing 
slope.  The  time  window  for  a  measurement  would  be  twice  the  time  from 
the  point  of  measurement  to  the  absorber  reflection.  One  could  determine 
the  point  of  measurement  by  shorting  the  G-TEM  cell  septum  to  the  shell 
where  the  sensor  would  be  placed  and  noting  where  the  reflection  occurs. 
A  less  exact  rule  of  thumb  is  to  divide  twice  the  distance  to  the  absorber 
by  the  propagation  velocity  in  air  of  1  ft/ns.  This  results  in  a  20-ns  win¬ 
dow  for  this  measurement  location. 

The  applied  waveform  is  formed  by  the  output  of  a  coaxial  charge  line 
connected  by  a  wetted  mercury  reed  switch  to  a  wave-shaping  circuit. 

The  wave-shaping  circuit  consists  of  parallel  resistors  before  and  after  a 
series  capacitor.  The  resistors  bleed  off  the  residual  charge  between 
pulses.  The  pulse  is  a  fast-rise  transient  with  an  exponential  fall  deter¬ 
mined  by  the  capacitor.  The  reed  switch  coil  is  operated  by  commercial 
60-Hz  ac  supplied  through  a  variac. 

The  pulser  output  was  applied  to  the  G-TEM  cell  as  a  source.  A  Tektronix 
7854  digital  sampling  oscilloscope  was  used  to  measure  the  signal.  A 
Tektronix  7M11  delay  line  was  needed  to  provide  the  necessary  75-ns 
signal  delay  after  triggering.  Spurious  triggering  was  occurring  on 
reflected  peaks  higher  than  the  incident  pulse  while  using  the  sensor 
signal  as  a  trigger  during  a  preliminary  measurement.  We  obtained  an 
accurate  trigger  by  splitting  the  pulser  output  and  using  half  for  a  trigger 
signal  and  half  for  the  G-TEM  source.  The  exact  configuration  for  apply¬ 
ing  the  signal  to  the  cell  was  measured  through  45  dB  attenuation  to 
protect  the  oscilloscope  input.  A  transfer  fimction  was  found  for  the 
attenuators,  oscilloscope  cables,  and  delay  line  used  in  the  measurement. 
The  pulser  output  tee  and  cable  to  the  G-TEM  cell  were  not  included  in 
the  transfer  function,  since  they  are  used  in  the  application  of  the  signal  to 
the  cell,  and  any  signal  deterioration  caused  by  these  components  should 
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Figxire  4.  TDR 
measurements  of 
G-TEM  cell  showing 
SO-Q  cable  open  and 
connect  to  cell- 


be  included  in  the  corrected  data.  The  applied  signals  were  varied  in 
amplitude  to  better  illustrate  the  different  sensor  responses  and  are 
plotted  in  figure  5.  The  pulses  are  all  close  to  4.5  ns  wide,  which  is  short 
enough  to  avoid  reflection  problems. 

Three  sensors  were  used  to  measure  each  electric  field  transient.  The 
sensors  were  identical,  except  for  having  lOx,  5x,  and  no  amplification. 
Each  sensor  was  placed  midway  between  the  septum  and  the  floor  of  the 
G-TEM  cell  at  the  0.9-m  septum  height.  The  sensor  cable  was  routed  in 
the  direction  of  field  propagation  for  a  meter  and  then  to  the  cell  wall 
feedthrough.  Care  was  taken  to  route  the  cable  perpendicular  to  the 
electric  field  when  possible. 

The  amplifiers  were  powered  through  the  center  conductor  by  a  dc  power 
supply  set  to  supply  60  mA  (around  5  V).  The  dc  voltage  was  blocked 
from  the  oscilloscope  input  by  a  bias  tee.  As  in  the  pulser  output  measure¬ 
ment,  the  sampling  oscilloscope  required  a  75-ns  delay  line  to  properly 
sample  the  signal.  All  the  components  between  the  oscilloscope  and  the 
sensor  were  characterized  in  a  transfer  ftmction.  We  made  a  transfer 
function  for  the  sensor  amplifier  using  two  bias  tees  and  the  pins  from  the 
sensor  with  no  amplifier,  mated  to  the  pins  of  the  sensor  with  an  ampli¬ 
fier.  The  transfer  functions  for  the  amplifiers  are  shown  in  figure  6. 

The  amplifiers,  although  specified  to  be  flat  to  4  GHz,  roll  off  at  aroimd 
2.5  GHz.  The  risetime  of  the  applied  pulses  is  averaged  to  be  135  ps, 
which  is  equivalent  to  an  upper  band  limit  of  2.6  GHz.  This  indicates  that 
some  amplifier  rolloff  may  affect  the  sensor  response.  A  correction  for  this 
rolloff  was  included  in  the  measurement  system  characteristics  correction. 
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There  were  measurements  made  with  each  sensor  dipole  inverted  to 
allow  elimination  of  common  mode  noise  coupled  to  the  sensor  cables. 
The  "down"  measurement  data  (dipole  ground  leg  down)  are  inverted 
and  added  to  the  "up"  data,  which  eliminate  the  common  mode.  The 
result  is  divided  by  two  to  compensate  for  the  signal  data  that  were 
added.  The  data  were  then  corrected  for  measurement  system  effects. 

Care  was  taken  to  only  invert  the  dipole  sensor  without  moving  the 
position  of  the  sensor  and  cable.  Care  was  also  taken  to  not  move  the 
sensor  and  cable  location  for  the  three  sensor  measurements  and  the  three 
field  levels  measured. 

Next,  we  calculated  the  electric  field  from  the  corrected  data  using  expres¬ 
sions  for  electrically  small  dipoles  [4].  An  electrically  small  dipole  has  an 
open  circuit  voltage  proportional  to  the  incident  electric  field  parallel  to 
the  anterma  in  the  time  or  frequency  domain, 

VJt)  =  2h,Ej,(t)  ,  (3) 


where  hg  is  the  equivalent  height  of  the  dipole.  A  relationship  for  the 
short-circuit  current  and  a  substitution  from  equation  (3)  results  in 


fsc(f)  ~ 


dt 


=  2h 


e'-a  f 


(4) 


where  is  the  antenna  capacitance.  In  relation  to  the  very  large  imped¬ 
ance  of  an  electrically  small  antenna,  the  50-Q  impedance  of  the  measure¬ 
ment  system  is  considered  to  be  a  short  circuit,  which  allows  the  use  of 
equation  (4)  to  relate  the  measurement  to  the  incident  field  [4].  A  factor  of 
50  is  introduced  from  V  =  I*R  relation  to  use  the  measured  voltage  data. 
Equation  (4)  is  integrated  on  both  sides  and  manipulated  to  solve  for  the 
electric  field. 


(5) 


where  Agq  is  called  the  equivalent  area  and  is  related  to  the  physical 
characteristics  that  define  the  dipole  and  antenna  capacitance.  The  cable 
impedance  Zg  is  50  Q  and  Ej,  is  the  electric  field  parallel  to  the  dipole.  The 
variables  defining  the  physical  characteristics  are  calculated  [5] 


C 

A  -  ‘’h 

^eq-  Eo'^e 


h  - 

^  r2-2  +  ln(4) 
2nhEQ 

Q-2-ln(4) 


(6) 


£2  =  2ta(f) 


5.  Results 


Figure  7.  Fields 
calculated  from 
measured  dipole 
voltage  compared 
with  245-V  G-TEM 
field  at  measurement 
location.  Comparison 
is  in  (a)  time  and 
(b)  frequency 
domains.  Constant 
factors  of  1.3  and  1.1 
are  used  to 
compensate  for 
inaccuracies  in  dipole 
model  predictions  for 
lOx  and  5x  sensors, 
respectively. 


where  h  is  the  half  dipole  length  and  a  is  the  dipole  radius.  The  results  of 
the  calculations  for  the  electric  field  of  the  lOx  and  the  5x  sensors  were 
multiplied  by  constants  of  1.3  and  1.1,  respectively.  This  was  to  match 
them  to  the  G-TEM  field  being  used  to  calibrate  the  sensors.  This  small 
adjustment  is  compensation  for  differences  between  ideal  dipole  models 
and  the  real  world.  The  MATLAB  routine  written  to  perform  this  series  of 
calculations  is  shown  in  appendix  C. 


The  measured  voltage  data  from  the  dipoles  placed  in  the  G-TEM  electric 
field  were  corrected,  integrated,  and  used  in  equation  (5)  to  calculate  the 
measured  electric  field.  Tlie  measured  electric  field  was  compared  with 
the  G-TEM  field  at  the  measurement  location.  The  results  are  shown  in 
figures  7,  8,  and  9. 
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Figure  8.  Fields 
calculated  from 
measured  dipole 
voltage  compared 
with  155-V  G-TEM 
field  at  measurement 
location.  Comparison 
is  in  (a)  time  and 
(b)  frequency 
domains.  Constant 
factors  of  1.3  and  1.1 
are  used  to 
compensate  for 
inaccuracies  in  dipole 
model  predictions  for 
lOx  and  5x  sensors, 
respectively. 


Frequency  (Hz) 


Figure  9.  Fields 
calculated  from 
measured  dipole 
voltage  compared 
with  60-V  G-TEM 
field  at  measurement 
location.  Comparison 
is  in  (a)  time  and 
(b)  frequency 
domains.  Constant 
factors  of  1.3  and  1.1 
are  used  to 
compensate  for 
inaccuracies  in  dipole 
model  predictions  for 
lOx  and  5x  sensors, 
respectively. 
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The  comparison  shows  an  improvement  in  the  high-frequency  character¬ 
istics  with  increased  amplification.  A  smaller  improvement  was  seen 
between  the  lOx  and  5x  sensors  because  they  both  increased  the  portion 
of  the  signal  lost  to  cable  attenuation  above  the  noise  floor.  A  much 
smaller  signal  would  be  needed  to  begin  degrading  the  amplified  sensor 
response  noticeably.  The  rolloff  of  the  amplifiers  at  2.5  GHz  causes  some 
loss  of  the  initial  E-dot  data,  because  the  upper  limit  of  the  differential 
field  frequency  content  is  twice  that  of  the  electric  field.  This  is  obvious 
because,  by  definition,  the  differential  field  has  its  peak  aroimd  the 
middle  of  the  risetime  of  the  electric  field.  As  mentioned  earlier,  the 
applied  pulses  have  frequency  content  up  to  2.6  GHz.  Using  the  reason¬ 
ing  above,  the  E-dot  signal  would  have  an  upper  frequency  content  of 
around  5  GHz.  The  frequency  content  between  2.5  and  5  GHz  was  not 
amplified  consistently,  and  the  transfer  function  for  the  amplifier  had  to 
be  applied  to  correct  inconsistencies.  The  cable  losses  present  at  the 
frequencies  not  amplified  decreased  some  of  the  signal  below  the  meas¬ 
urement  sensitivity  of  the  equipment.  This  lost  information  would  cause 
the  applied  transfer  function  to  lose  accuracy.  This  would  accoimt  for  the 
inability  to  measure  the  same  peak  value  as  the  G-TEM  field  and  the  loss 
of  some  of  the  structure  in  the  transient.  Even  with  the  amplifier's  capa¬ 
bilities  being  exceeded,  the  data  show  an  increase  in  accuracy  over  the 
nonamplified  dipole. 

Figure  10  is  an  E-dot  data  plot  of  the  60  V/m  peak  field  measiured  with 
the  5x  dipole.  The  upper  plot  of  figure  10  also  shows  the  common  mode 
noise  caused  by  fields  coupling  to  the  sensor  cable.  The  common  mode 
signal  is  present  at  the  dipole,  where  it  combines  with  the  ground  leg 
portion  of  the  dipole  signal.  Inverting  the  dipole  inverts  the  dipole  re¬ 
sponse,  while  the  common  mode  noise  remains  the  same.  This  allows  us 
to  remove  the  noise  by  inverting  the  mverted  signal  digitally,  adding,  and 
dividing  by  two  to  correct  the  amplitude.  The  result  of  this  process  is 
shown  in  the  lower  plot. 

Figure  11  compares  the  E-dot  data  from  the  lOx,  5x,  and  lx  sensors.  With 
the  extrapolation  of  the  nonamplified  sensor  response  peaks  from  field 
transients  with  peaks  of  245  and  155  V/m  in  the  bottom  plot  of  figure  11, 
the  measurement  of  the  field  transient  with  a  peak  of  60  V/m  should 
provide  a  sensor  response  peak  of  aroimd  8  mV.  A  measurement  of  the 
60  V/m  peak  field  transient  with  the  unamplified  dipole  was  not  possible 
with  the  available  scope  sensitivity. 

Assuming  a  minimum  usable  transient  peak  amplitude  of  10  mV,  we  can 
use  the  E-dot  sensor  response  peaks  in  the  top  and  middle  plots  of  figure 
11  to  determine  the  field  level  that  would  result  in  a  10-mV  peak.  If  we 
assume  linearity,  the  ratio  of  a  measured  E-dot  sensor  peak  response  to 
the  G-TEM  field  times  the  minimum  measurable  peak  response  of  10  mV 
should  provide  the  minimum  field  measurable  by  that  sensor.  The  appli¬ 
cation  of  this  proportionality  to  the  5x  sensor  data  results  in  a  minimum 
measurable  field  peak  of  13  V/m  from  an  average  of  the  three  measure¬ 
ments  with  the  sensor.  The  lOx  sensor  minimum  measurable  field  peak  is 
8  V/m. 
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Figiire  10.  Sensor 
response:  (a)  Raw 
dipole  response  for 
ground  leg  down  and 
up.  Measurement  was 
made  with  5x  sensor 
placed  in  60  V/m  peak 
electric  field. 

(b)  Result  of 
combining  two  data 
sets  to  remove  noise. 


Figure  11.  Results  of 
E-dot  data  with 
common  mode  noise 
removed  in  G-TEM 
fields  with  (a)  lOx, 
(b)  5x,  and  (c)  no 
amplifier.  An 
individual  sensor's 
response  to  different 
G-TEM  fields  is 
compared  in  each 
plot. 
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The  peak  amplitudes  of  the  lOx  and  5x  sensors  would  differ  by  a  factor  of 
two  if  the  amplification  was  flat  in  the  frequency  domain.  We  find  a  factor 
of  aroimd  1.7  by  inspecting  for  the  difference  in  peak  amplitudes  from  the 
lOx  and  5x  sensors.  This  agrees  with  the  value  of  1.75  we  find  by  averag¬ 
ing  the  result  of  dividing  the  transfer  function  for  the  lOx  amplifier  by  the 
transfer  function  for  the  5x  amplifier  and  restricting  the  data  to  below 
3  GHz.  The  lOx  sensor  data  have  a  peak  amplitude  factor  of  arormd  8 
with  respect  to  the  nonamplified  sensor  response. 

As  the  measured  amplitudes  near  the  lower  bound  of  sensitivity,  there  is 
loss  of  structure  in  the  time-domain  signal  and  a  corresponding  loss  of 
frequency  content.  The  structure  loss  occurs  when  the  lower  amplitude 
peaks  are  lost  in  the  noise  of  a  time-domain  measurement.  If  the  amplifi¬ 
cation  at  the  sensor  is  constant  in  the  frequency  domain,  the  minimum 
measurable  unamplified  signal  is  increased  by  the  amplification  factor. 
The  amplified  signal  will  still  suffer  degradation  as  it  nears  the  improved 
lower  botmd  of  sensitivity. 

A  resonance  in  the  frequency  domain  can  be  seen  in  the  frequency  re¬ 
sponse  plots  of  the  sensors  (fig.  6-8)  at  around  2.1  GHz.  This  is  a  result  of 
the  construction  of  the  sensors  and  can  be  minimized  with  modifications. 
The  suspect  area  would  be  the  dipole  gap  where  stray  capacitances  or 
inductances  could  occur.  The  fact  that  the  resonance  is  not  as  dramatic  for 
the  measurements  in  the  lowest  field  levels  could  indicate  a  higher  noise 
floor. 

Fabrication  methods  could  be  improved  for  the  amplifier  installation. 

Heat  limits  were  probably  exceeded  when  the  coax  shield  was  soldered. 
This  could  explain  the  loss  of  amplification  above  2.6  GHz  when  the 
specifications  are  4  GHz.  Same-size  amplifiers  with  upper  limits  at  6  GHz 
are  manufactured  by  the  same  manufacturer.  A  stable  65-mA  current 
controlled  voltage  source  could  be  made  to  reduce  small  errors  due  to 
inconsistent  voltage  supply. 

The  sensor  response  in  the  frequency  domain  with  the  measurement 
system  and  amplification  removed  is  shown  in  figure  12.  The  rolloff  at  the 
high-frequency  end  is  the  result  of  filtering  to  remove  large  oscillations 
caused  by  noise. 


6.  Conclusions 

The  amplifier  incorporated  into  the  sensor  improves  the  measurement 
system  sensitivity  and  accuracy  by  the  performance  of  the  amplifier.  We 
avoid  any  signal  degradation  due  to  measurement  system  characteristics 
by  amplifying  the  signal  before  it  can  be  affected.  An  improvement  in  the 
fabrication  methods  and  choice  of  amplifier  that  produces  flat  response  in 
the  frequency  domain  being  measured  can  improve  the  accuracy  of  the 
sensor.  The  sensors  used  in  the  experiment  show  improvement  in  sensi¬ 
tivity  and  accuracy  over  the  nonamplified  sensor,  despite  inconsistent 
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Figure  12.  Transfer 
functions  for  three 
sensors  used  in 
experiment,  corrected 
for  frequency 
characteristics  of 
measurement  system 
and  amplifiers  in 
sensor  with  (a)  lOx, 
(b)  5x,  and  (c)  no 
amplifier.  Each  plot 
contains  result  of 
indicated  sensor's 
response  to  G-TEM 
field  convolved  with 
the  G-TEM  field  in 
frequency  domain. 


Frequency  (GHz) 


Frequency  (GHz) 

amplification  in  the  measurement  frequency  band.  An  increase  in  sensi¬ 
tivity  was  shown,  from  above  60  to  13  V/m  for  the  5x  sensor  and  8  V/m 
for  the  lOx  sensor. 

The  feature  of  pin  connectors  in  the  sensor  allows  for  flexibility  in  its  use. 
For  example,  an  H-field  sensor  could  be  made  using  a  small  loop  in  place 
of  the  dipole.  For  cw  measurements,  we  could  enhance  the  sensitivity  of 
the  sensor  by  sizing  the  antenna  to  be  resonant.  This  would  vastly  im¬ 
prove  the  antenna  response  above  that  of  the  electrically  small  antenna 
necessary  for  time  domain  measurements. 

The  size  of  the  sensor  is  minimized  to  the  dipole  and  the  semirigid  coaxial 
measurement  cable.  This  minimum  size  reduces  the  perturbation  of  the 
field  being  measured  by  the  presence  of  the  sensor.  The  s)anmetry  of  the 
sensor  allows  scattered  fields  incident  from  above  and  below  the  sensor 
to  be  measured  without  interference.  Ground  plane  sensors  cannot 
accurately  perform  these  measurements.  The  symmetry  also  permits  the 
sensor  to  be  inverted  to  provide  the  data  that  allow  removal  of  common¬ 
mode  noise  from  the  sensor  response. 
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More  complicated  circuits  could  be  made  small  enough  to  remain  in  the 
electrically  small  category  and  provide  a  high  impedance  for  direct 
electric  field  measurements,  antenna  loading,  etc.  The  miniaturization  of 
electronic  components  makes  what  was  improbable  a  few  years  ago  a 
possibility  today. 
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Appendix  A. — ^MATLAB  Function  CALSTEP 

This  appendix  gives  the  MATLAB  routine  M-file  calstep.m,  which  is  the 
source  for  the  function  calstep.  The  function  calstep  uses  measured  time- 
domain  data  for  component  input  and  output  to  create  a  component 
transfer  function  by  emplo)dng  an  optimal  deconvolution  method  in  an 
iterative  fashion.  Progress  toward  a  solution  is  determined  using  a 
graphic  interface  and  knowledge  of  the  frequency  content  of  the  applied 
waveform. 

%  b  =  the  input  signal 
%  a  =  the  output  signal 
%  c  =  the  transfer  function 
% 

%  fimction  ccc  =  calstep  (b,  a) 

% 

%  check  the  inputs  to  see  if  they  are  a  power  of  two  and 

%  equally  spaced.  If  not,  interpolate  to  a  power  of  two 

%  which  also  matches  the  waveforms  point  for  point. 

% 

function  ccc  =  calstep (b, a) 

%  First,  ensure  the  first  point  is  zero  time  and  amplitude 
if  a(l,l)  ~=  0,  a(l,l)  =  0.0;  end 

if  a(l,2)  ~=  0,  a(l,2)  =  0.0;  end 

if  b(l,l)  -=  0,  b(l,l)  =  0.0;  end 

if  b(l,2)  -=  0,  b(l,2)  =  0.0;  end 

nb  =  length (b) ; 
na  =  length (a) ; 
tb  =  b(nb, 1) ; 
ta  =  a (na, 1) ; 
intpl=0 ; 

xx=diff (b( : , 1) ) ; 

if  (1  -  min(xx)  /  max(xx))  >  .01 

disp('  Need  equally  spaced  points  in  the  input  data') 
disp('  The  points  will  be  interpolated  to  the  smallest  dt') 
intpl=l ; 
end 

dtmin  =  min(xx); 
xx=diff (a( : , 1) )  ; 

if  (1  -  min(xx)  /  max(xx))  >  .01 

disp('  Need  equally  spaced  points  in  the  output  data') 
disp('  The  points  will  be  interpolated  to  the  smallest  dt') 
intpl=l; 
end 

dta  =  min(xx); 
if  dta  <  dtmin 
dtmin  =  dta; 
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end 

ninax=nb  ; 

while  log2  (nitiax)  >  fix(log2  (nmax)  ) 

disp('Need  a  power  of  two  data  points,  will  interpolate') 

disp(  [  'Total  points  now  are  '  ,nuin2str  (nmax)  ]  ) 

nmax= input ( 'How  many  points  do  you  want?') 

dtmin  =  tb/ (nmax-1) ; 

intpl=l; 

end 

tmax=tb; 

t  =  dtmin  *  (0 :nmax-l) ' ; 
if  tb  >  ta 

disp('  Need  equal  length  data  sets,  input  is  longer') 
disp('  The  data  will  be  interpolated  to  the  shortest  length') 
nmax  =  fix(ta/dtmin)  +  1; 
if  log2 (nmax)  >  fix (log2 (nmax) ) 

disp('Need  a  power  of  two  data  points,  will  interpolate') 
nmax  =pow2 ( f ix ( log2 (nmax) ) +1 ) ; 
dtmin  =  ta/ (nmax-1) ; 

end 

tmax  =  (nmax-1)  *  dtmin; 
t  =  dtmin  *  ( 0 : nmax- 1 ) ' ; 
intpl=l ; 

end 

if  ta  >  tb 

disp('  Need  equal  length  data  sets,  output  is  longer') 
disp('  The  data  will  be  interpolated  to  the  shortest  length') 
nmax  =  f ix ( tb/dtmin) +1 ; 
if  log2 (nmax)  >  f ix (log2 (nmax) ) 

disp('Need  a  power  of  two  data  points,  will  interpolate') 
nmax  =pow2 (fix (log2 (nmax) ) +1) ; 
dtmin  =  tb/ (nmax-1) ; 

end 

tmax  =  (nmax-1)  *  dtmin; 
t  =  dtmin  *  (0:nmax-l)'; 
intpl=l ; 

end 

if  intpl  ==  1 

sa  =  interpl (a( : , 1) ,a( : ,2) ,t) ; 
sb  =  interpl (b ( : , 1 ) , b ( : , 2 ) , t ) ; 
else 

sa  =  a ( : , 2 ) ; 
sb  =  b ( : , 2 ) ; 

end 

subplot (2,1,1) 

plot (a ( : ,1) ,a( : ,2) , t,sa) 

title ( 'Output ' ) 

subplot (2, 1,2) 

plot (b( : , 1) ,b( : , 2 ) , t, sb) 
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title ( 'Input' ) 
figure (1) 
pause (2 ) 

% 

%  The  Nicolson  ramp  method  and  the  Gans-Nahman  inverse  waveform 
%  method  are  used  here. 

%  The  Nicolson  ramp  method  uses  a  ramp  "ra"  created  using  the  number 
%  of  points  "na"  and  the  amplitude  of  the  final  point: 
n  =  0:l:nmax-l; 

ra  =  (n/(nmax-l)  *  sa(nmax))'; 

rb  =  (n/(nmax-l)  *  sb(nmax))'; 

sra  =  sa  -  ra; 

srb  =  sb  -  rb; 

subplot (2,1,1) 

plot (t, sra) 

title ( 'Output  ramped  to  O') 
siibplot  (2,1,2) 
plot ( t, srb) 

title ( ' Input  ramped  to  O') 
figure (1) 
pause (2 ) 

sram  =  fft(sra); 
srbm  =  fft(srb); 

f  =  (0:l/t(nmax) :l/((t(2)-t(l))*2))'; 
subplot (2,1,1) 
sam=abs (f f t (sa) ) ; 

loglog ( f , abs ( sram ( 1 : length ( f ) ) ) , f , sam ( 1 : length ( f ) ) ) 
title ('ramp  FFT  vs  FFT  for  output') 
subplot (2,1,2) 
sbm=abs (f ft (sb) ) ; 

loglog ( f , abs ( srbm ( 1 : length ( f ) ) ) , f , sbm ( 1 : length ( f ) ) ) 
title('ramp  FFT  vs  FFT  for  input') 
figure (1) 
pause (2 ) 

% 

%  The  Gans-Nahman  inverse  waveform  method  appends  an  inverse 
%  waveform  "u"  to  the  last  point  of  the  original  waveform. 

% 

ua  =  sa(nmax)  -  sa; 
ub  =  sb(nmax)  -  sb; 
sua  =  [sa  ;  ua] ; 
s\ib  =  [sb  ;  ub]  ; 

tu  =  0:(2  *  tmax/(2  *  nmax  -  1)):2  *  tmax; 
subplot (2,1,1) 
plot ( tu, sua) 

title ( 'Output  for  Gans-Nahman  invers  waveform  FFT') 
subplot (2, 1,2) 
plot ( tu, sub) 

title ('Input  for  Gans-Nahman  invers  waveform  FFT') 
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figure (1) 
pause (2 ) 

suam  =  fft{sua)  *  0.5; 
subrti  =  f  ft  (sub)  *  0.5; 

f  =  (0:l/tu{2  *  ninax)  :!/ (  {tu(2) -tu(l)  )  *2)  )  '  ; 

% 

%  Replace  the  even  function  series  of  the  Gans-Nahman  method 
%  with  the  values  calculated  in  the  Nicolson  method.  This 
%  would  replace  the  erronious  even  values  caused  by  the 
%  frequency  spacing.  The  DC  value  is  correct  for  the  Gans-Nahman 
%  method  and  will  remain. 

% 

K=3; 

cma=suam; 

for  k  =  l:nmax-l 

cma{K)  =  sram(k+l) ; 

K=K+2 ; 

end 

K=3; 

cmb=subm; 

for  k  =  l:nmax-l 

cmb(K)  =  srbm{k+l) ; 

K=K+2 ; 

end 

siibplot  (2,1,1) 

loglog(f ,abs (cma(l:nmax) ) ) 

title ( 'Combined  Gans-Nahman  inverse  and  ramp  method  FFT  for  Output') 

subplot (2,1,2) 

loglog (f , abs (cmb (1 :nmax) ) ) 

title ( 'Combined  Gans-Nahman  inverse  and  ramp  method  FFT  for  Input') 
figure (1) 
pause (2 ) 

% 

%  Perform  the  division  of  output  over  input  to  obtain  the  transfer 

%  function  for  the  common  period  and  the  smallest  dt  of  the  input 

%  and  output  functions . 

% 

%  This  is  done  using  a  deconvolution  algorithm. 

% 

%  or  by  simply  dividing 
% 

c  =  cma  . /  cmb ; 
c(l)  =  0; 

%  assiimes  a  DC  content  of  zero  in  the  transfer  function 
% 

siibplot  (1,1,1)  ; 

rta=risetime (a) ; 
pause (2 ) 
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rtb=risetiine  (b)  ; 
pause ( 2 ) 
rtfa= . 35/rta; 
rtfb= . 35/rtb; 
close (1) 

nrtfa=max(find{f  <  rtfa))+l; 
nrtfb=]nax(find(f  <  rtfb))+l; 

crta=  [inin(abs  (c  (2  :nrtfa)  )  )  /lO  10*max(abs  (c  (2  rnrtfa)  )  )  ]  ; 
crtb=  [itiin{abs  (c  (2  :nrtfb)  )  )  /lO  10*inax{abs  (c  (2  :nrtfb)  )  )  ]  ; 

%  using  the  deconvolution  algorithm 
%  With  an  iteration  program  to  specify  G 
f  =  (0:l/tu(2  *  nmax) :1/ ( (tu(2)-tu(l) ) ) ) ' ; 

maxfreq=l/ (dtmin*10) ; 
fnum=min(find(f  >=  maxfreq) ) ; 

G  =  input ('Enter  a  value  for  the  iteration  '); 

% 

figure (1) 
subplot (1,1,1) ; 
while  G  ~=  0 

CFILTER=conj  (cmb)  ./  (abs(cinb)  2  +  G)  ; 

CF  =  cma  .*  CFILTER; 

CF(1)  =  0; 

%  assximes  a  DC  content  of  zero  in  the  transfer  function 

loglog (f (1 :nmax) , abs (CF (1 :nmax) ) , f (1 :nmax) , abs (c (1 :nmax) ) , . . . 
[rtfa  rtfa] , crta, , [rtfb  rtfb] , crtb, ' : ' ) 
axis([f(l)  f (nmax)  crta] ) ; 
pause (1) 

title ([ 'Enter  value  or  0  to  quit.  Base  Risetime  freq  = 

'  , num2str (rtfb) ] ) 

GG=G; 
figure (1) 

G  =  input ('Enter  a  value  for  the  iteration,  0  to  quit  ' ) ; 
end 


% 

% 

loglog (f (1 :nmax) , abs (CF (1 :nmax) ) , f (1 : nmax) , abs (c (1 :nmax) ) ) 

title ([ 'Convolution  with  constant  '  ,n\im2str (G)  ,  '  vs  constant  =  0']) 

pause (2 ) 

loglog  (f  (1 :  fnum)  ,  abs  (CF  (1 :  fn-om)  )  ,  [rtfa  rtfa]  ,  crta,  '  — '  ,  [rtfb 
rtfb] , crtb, ' :  '  ) 

title ('Cal  curve  plotted  to  the  frequencies  with  >10  points  in  the 
risetime' ) 
ccc= [f  CF] ; 
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This  appendix  gives  the  MATLAB  routine  M-file  cal.m,  which  is  the 
source  for  the  function  cal.  The  function  cal  allows  measured  data  to  be 
"calibrated"  or  to  have  the  effects  of  the  frequency-domain  characteristics 
of  the  measxurement  system  components  contained  in  a  transfer  function 
to  be  removed.  The  transfer  function  is  created  with  calstep.m,  the 
MATLAB  fimction  given  in  appendix  A  of  this  report. 


%  This  function  takes  a  transfer  fimction  "tfn"  that  has  been  ob¬ 
tained  using 

%  the  function  calstep  on  two  equal  length  time  files  that  describe 
the 

%  response  of  the  system  with  and  without  the  portion  being  cali¬ 
brated  and 

%  applies  it  to  the  measured  signal  "s".  "s"  must  be  the  same  length 
as  the 

%  transfer  function  time  files 
% 

% 

%  function  sysresp  =  cal (tfn, s) 

% 

% 

%  check  the  inputs  to  see  if  they  are  a  power  of  two  and  equally 
%  spaced.  If  not,  interpolate  to  a  power  of  two  which  also 

%  matches  the  waveforms  point  for  point. 

% 

function  sysresp  =  cal ( tfn, s) 

% 

%  First,  ensure  the  first  point  is  zero  time  and  amplitude 
if  s(l,l)  ~=  0,  s(l,l)  =  0.0;  end 
if  s(l,2)  ~=  0,  s(l,2)  =  0.0;  end 
sa=s ( : , 2 ) ; 
t=s ( : , 1)  ; 
ns  =  length(sa); 

% 

% 

xx=di f  f ( s ( : , 1 )  )  ; 

if  (1  -  min{xx)  /  max(xx))  >  .01 

disp('  Need  equally  spaced  points  in  the  input  data') 
disp('  The  points  will  be  interpolated  to  the  smallest  dt') 
intpl=l; 
pause (2 ) 
end 

dtmin  =  mean(xx) ; 
if  log2(ns)  >  fix(log2 (ns) ) 

disp('Need  a  power  of  two  data  points,  will  interpolate') 
ns  =pow2 (fix(log2 (ns)  ) -1)  ; 
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dtmin  =  s (length (s) , 1) / (ns-1) ; 

t  =  dtmin  *  (0:ns-l)'; 

sa  =  interpl { s ( : , 1 ) , s ( : , 2 ) , t ) ; 

end 


%  The  Nicolson  ramp  method  and  the  Gans-Nahman  inverse  waveform 
%  method  are  used  here. 

%  The  Nicolson  ramp  method  uses  a  ramp  "ra"  created  using  the  n\imber 
%  of  points  "n"  and  the  amplitude  of  the  final  point: 

% 

n  =  0:l:ns-l; 

ra  =  (n/(ns-l)  *  sa(ns))'; 
sra  =  sa  -  ra; 

plot (s(:,l),s(:,2),t, sra , ' . ' ) 

title ( 'Original  data  versus  data  ramped  to  zero') 

figured) 

pause (2 ) 

sram  =  fft(sra); 

f  =  (0:l/t(ns) :l/((t(2)-t(l))*2))'; 
sam=abs (f f t (sa) ) ; 

loglog ( f , abs ( sram ( 1 ; length ( f ) ) ) , f / sam ( 1 : length ( f ) ) ) 
title ( 'Nicolson  ramp  FFT  vs  regular  FFT') 
figure (1) 
pause (2 ) 

% 

%  The  Gans-Nahman  inverse  waveform  method  appends  an  inverse 
%  waveform  "u"  to  the  last  point  of  the  original  waveform. 

% 

ua  =  sa(ns)  -  sa; 
sua  =  [sa  ;  ua] ; 
tu  =  dtmin  *  (0:2*ns-l); 
plot (tu, sua) 

title ('The  Gans-Nahman  inverse  waveform  to  be  FFTed' ) 
figure (1) 
pause (2 ) 

suam  =  f  f  t ( sua )  *  0.5; 

f  =  (0:l/tu(2  *  ns) :1/ ( (tu(2)-tu(l) ) *2) ) ' ; 


% 

%  Replace  the  even  function  series  of  the  Gans-Nahman  method  with 
the 

%  values  calculated  in  the  Nicolson  method.  This  would  replace  the 
%  erronious  even  values  caused  by  the  frequency  spacing.  The  DC 
value 

%  is  correct  for  the  Gans-Nahman  method  and  will  remain. 

% 

K=3; 

sf=suam; 

for  k  =  l:ns-l 
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sf(K)  =  sram(k+l) ; 

K=K+2 ; 

end 

subplot (1,1,1) 

loglog ( f , abs ( sf ( 1 : length ( f ) ) ) , tf n ( : , 1 ) , abs ( tf n { : , 2 ) ) ) 

title ( 'Resulting  combination  of  signal  FFTs  vs  transfer  function') 

figure (1) 

pause (2 ) 

% 

% 

% 

tf=tfn( : ,2) ; 

% 

%  Perform  the  division  of  output  over  input  to  obtain  the  transfer 

%  function  for  the  common  period  and  the  smallest  dt  of  the  input 

%  and  output  functions . 

% 

%  This  is  done  using  a  deconvolution  algorithm. 

% 

%  or  by  simply  dividing 
% 

c  =  sf  ./  tf; 
rt=risetime (s) ; 
pause (2) 
rtf= . 35/rt ; 
close (1) 

nrtf=max(find(f  <  rtf))+l; 

crt= [min (abs (c (2 :nrtf ) ) ) /lO  10*max(abs (c (2 :nrtf ) ) ) ] ; 

% 

%  using  the  deconvolution  algorithm 
%  With  an  iteration  program  to  specify  G 

maxfreq=l/(dtmin*10); 
fnum=min(find(f  >=  maxfreg) ) ; 

G  =  input ('Enter  a  value  for  the  iteration'); 

% 

subplot (1,1,1)  ; 
while  G  ~=  0 

CFILTER=conj  (tf )  ./  (abs(tf)  . 2  +  G)  ; 

CF  =  sf  .*  CFILTER; 

loglog(f (l:ns) , abs (CF (1 :ns) ) ,f (2:ns) ,abs(c(2:ns) ) , [rtf 
rtf] , crt , ' : ' ) 

axis([f(l)  f(ns)  crt] ) ; 
pause (1) 

title ([ 'Enter  value  or  0  to  quit.  Risetime  freq  = 

'  ,nvun2str  (rtf )  ]  ) 

GG=G; 
figure (1) 
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G  =  input ('Enter  another  value  for  the  iteration,  0  to  quit  ' )  ; 
end 

loglog (f (1 : fnum) , abs (CF (1 : fnum) ) , [rtf  rtf] , crt , ' : ' ) 
sysres  =  ifft(CF); 

title { 'Risetime  Frequency  is  dotted  line') 
sysresp= [t  real (sysres (1 : length(t) ) ) +ra] ; 
figure (1) 
pause 

disp  (['The  maximum  imaginary  part  of  the  IFFT  is:', 
num2str (max ( imag ( sysres ) ) ) ] ) 
plot ( sysresp ( : , 1 ) , sysresp ( : , 2 ) ) 
figure (1) 
pause (2 ) 
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Appendix  C. — Data  Correction  Algorithm 

This  is  the  algorithm  used  to  correct  the  data  for  measurement  system 
characteristics,  remove  common-mode  noise,  and  calculate  the  field  from 
the  dipole  dimensions.  The  other  routines  called  in  this  file  are  included 
at  the  end  of  this  appendix. 


%  This  routine  accelerates  the  process  of  using  the  dipole 
%  inverted  and  not  inverted  signal  to  remove  common  mode  cable 
%  signals.  The  measurement  must  be  calibrated  to  remove 
%  measurement  system  characteristics.  The  antenna  amplifier  at 
%  the  dipole  can  be  treated  as  a  constant  or  the  characteristcs 
%  can  be  included  in  with  the  measurement  system.  The  signal  with 
%  the  ground  leg  up  is  shifted  to  remove  DC  of f set, inverted, 

%  shifted  to  match  risetimes,  added  together  and  divided  by  2, 

%  corrected,  integrated,  then  the  E-field  proportional  to  the 
%  measured  voltage  is  calculated  from  dipole  dimensions  for  an 
%  electrically  small  antenna. 

% 

%  input  needed  is: 

% 

%  sup  =  measured  input  signal  with  dipole  shield  leg  up 

%  sdn  =  measured  input  signal  with  dipole  shield  leg  down 

%  tcbl  =  transfer  function  of  measurement  system  previously 

determined 

%  rdipole  =  dipole  radius  in  meters 

%  hdipole  =  dipole  half  height  in  meters 

% 

%  function  aaa=cdipole (sup,  sdn,  tcbl,  rdipole,  hdipole); 

% 

% 

function  aaa=cdipole ( 1 ,  sdn,  tcbl,  rdipole,  hdipole,  amp) ; 

%  Remove  cable  losses  using  transfer  functions  generated  with 
calstep 

%  First  remove  any  DCoffset 
sup=dcof f set (sup) ; 

sdn=dcof f set (sdn) ; 

%  invert  the  reversed  measurement 

aaa=sup; 

aaa( : ,2)=-sup( : ,2) ; 

disp('Now  you  must  line  up  the  rise  times  of  the  signals') 
pause (2 ) 

aaa=shif ty (aaa, sdn) ; 


27 


Appendix  C 


%  add  the  data  and  divide  by  two,  this  removes  common  mode 
noise  bbb=aaa ; 

bbb( : , 2 ) = (aaa ( : , 2 ) +sdn ( : , 2 ) ) /2 ; 

%  set  first  and  last  point  to  zero 

bbb(l,2)=0; 

bbb(l,2)=0; 

ccc=cal (tcbl,  bbb); 

ddd=shif tint (ccc) ; 
pause (2 ) 

%  now  the  antenna  Aeq  is  calculated 
close (1) 
paused) 

eps0=8 . 854e-12 ; 

omega=2  * log ( 2  *hdipole/rdipole ) ; 

disp([  'Antenna  fatness  factor  is  nijim2str  (omega)  ] ) 

heq=hdipole* (omega-1) / (omega-2 +log (4) ) ; 

disp([  'Antenna  equivalent  height  is  ',  num2str (heq) ] ) 

Ca=2*pi*hdipole*eps0/ (omega-2 -log (4) ) ; 

disp([  'Antenna  capacitance  is  ',  num2str (Ca) ] ) 

Aeq=Ca / eps  0  *  heq ; 

disp([  'Antenna  equivalent  area  is  ',  num2str (Aeq) ] ) 
pause (3 ) 

%  calculate  incident  field 
eee=ddd; 

eee ( : , 2 ) =ddd( : ,2) / (eps0*Aeq*50) ; 
plot  (eee  {•.,!),  eee  ( :  ,  2 )  ) 

title ([ 'Measured  E-field  with  dipole  ']) 
xlabel ( ' Seconds ' ) 
ylabel ( 'Volts/Meter' ) 

%***  DCOFFSET.M  *** 

% 

%  This  function  shifts  the  Y  values  up  or  down  by  the  mean  of 
%  the  baseline  data.  The  amount  of  baseline  used  is  determined 
%  visually. 
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%  The  data  must  be  two  columns  of  X,  Y  values 
% 

%  b  =  the  curve  to  be  shifted 
%  a  =  the  baseline  to  be  averaged 
% 

function  sh  =  dcoffset(b) 
a=b; 

plot (b(:,l),b(:,2),a(:,l),a(:,2)) 

title ('Input  the  begining  and  end  points  in  the  baseline  to  be  used 
') 

grid 

figure (1) 

disp( '*******  Enter  end  =  1  to  quit  ★******') 

ns  =  input ('Input  the  begining  point  in  the  baseline  to  be  used  '); 
ne  =  input ('Input  the  end  point  in  the  baseline  to  be  used  ' ) ; 

while  ne  >  1 

if  ns  ==  0,  ns  =  1;  end 
if  ne  ==  0,  ne  =  1;  end 
a  =  b(ns:ne,:); 

plot(b(:,l) ,b(:,2),a(:,l),a(:,2)) 

title ('Input  the  begining  and  end  points  in  the  baseline  to  be 
used' ) 

grid 

figure (1) 
nns=ns ; 

ns  =  input ('Input  the  begining  point  in  the  baseline  to  be  used 

ne  =  input ('Input  the  end  point  in  the  baseline  to  be  used  '); 
end 

ave=mean ( a ( :  ,  2 ) )  ; 
sh=b  ; 

sh ( : , 2 ) =b ( : , 2 ) -ave ; 

sh(l :nns , 2 ) =sh(l :nns ,2)*0; 

plot ( sh ( : , 1 ) , sh ( : , 2 ) ) 

title ( 'Shifted  to  eliminate  DC  offset') 
figured) 

%***  SHIFTINT.M  *** 

%  this  routine  combines  shiftx.m  and  funint.m  into  a  loop  to  allow 
the 

%  interactive  removal  of  DC  bias 
% 

% 

%  This  function  shifts  the  Y  values  up  or  down  while  plotting  the 
%  reference  and  shifted  curves  for  comparison 
%  The  data  must  be  two  col\imns  of  X,  Y  values 
% 
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%  b  =  the  curve  to  be  shifted 
% 

function  c  =  shiftint(b) 
a=b; 

c=funint (b) ; 
plot (c ( :,l),c(:,2)) 

title('The  integrated  unshifted  curve') 
figure (1) 
pause 

title (' ***Enter  value  to  shift  curve  +up  or  -down***') 
pause (2 ) 

plot(b( :,l),b(:,2),a(:,l),a(:,2)) 

title ('Enter  value  to  shift  curve  +up  or  -down'); 

figure (1) 

grid 

ns  =  input ('Input  the  value  to  shift  y  values  up  or  down  '); 
sh  =  b; 
while  ns  ~=  0 

sh ( : , 2 )  =  ns  +  b ( : , 2 ) ; 
c=funint (sh) ; 
plot(c( ; ,1) ,c( : ,2) ) 

title ('The  integrated  shifted  curve.  Enter  value  to  shift  again') 
figure (1) 

ns  =  input ('Input  the  value  to  shift  y  values  up  or  down  ' ) ; 
end 

c=funint (sh) ; 
plot(c( : ,1) ,c(: ,2) ) 

title('The  final  integrated  shifted  curve') 
figure (1) 

%***  SHIFTY. M  *** 


% 

%  This  function  shifts  the  Y  values  to  the  right  while  plotting  the 
%  reference  and  shifted  curves  for  comparison 
%  The  data  must  be  two  columns  of  X,  Y  values 
% 

%  b  =  the  curve  to  be  shifted 
%  a  =  the  reference  curve 
% 

function  sh  =  shifty (b,  a) 
plot (b(:,l),b(:,2),a(:,l),a(:,2)) 

title( 'Enter  #  of  spaces  to  shift  yellow  curve  +right  or  -left'); 
figure (1) 

nstog  =  input ('Enter  #  of  spaces  to  shift  yellow  curve  +right  or  - 
left'  )  ; 

np  =  length (b) ; 
while  nstog  ~=  0 
ns=nstog; 
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sh(:,l)  =  b(:,l)+(b(2,l)*(ns-l)); 

Sh(: ,2)  =  b(: ,2) ; 

plot ( sh ( : , 1 ) , sh ( : , 2 ) , a ( : , 1 ) , a ( : , 2 ) ) 

title('#  of  spaces  to  shift  yellow  curve  +right  or  -left'); 
figure (1) 

nstog  =  input ('Enter  #  of  spaces  to  shift  yellow  curve  +right  or  - 
left' ) ; 
end 

sh{:,l)=b(:,l); 
sh(:,2)  =  b(:,2)*0; 
if  ns  >  0 

sh(ns+l :np, 2 )  =  b ( 1 : (np-ns ) , 2 ) ; 
else 

ns=abs (ns) ; 

sh(l : (np-ns) , 2 )  =  b(ns+l;np,2); 

sh(np-ns+l :np, 2)  =  sh (np-ns+1 :np, 2 ) *0+sh (np-ns , 2 ) ; 

end 

plot ( sh ( : , 1 ) , sh ( : , 2 ) , a ( : , 1 ) , a ( : , 2 ) ) 
title( 'Final  shifted  curve') 
figure (1) 

%***  RISETIME. M  *** 

%  This  routine  finds  the  peak  and  computes  the  10  to  90  percent 
risetime 
% 

%  Input  the  signal  variable  -  s 
%  function  a=risetime (s) 

function  dxr=risetime ( s ) 

% 

plot (s ( :,l),s(:,2)) 

title ('Place  cursor  at  begining  of  risetime  peak,  click  left  button') 
aaa=ginput ( 1 ) ; 
xs=aaa(l)  ; 

title ('Place  cursor  at  end  of  risetime  peak,  click  left  button') 
aaa=ginput ( 1 ) ; 
xe=aaa (1); 

ns=max(find(s ( : , 1)  <=  xs)); 

ne=max(find(s ( : , 1)  <=  xe) ) ; 

plot (s ( : , 1) , s ( : , 2 ) , s (ns :ne, 1) , s (ns  ;ne,  2 ) ) 

title (' Selected  area  in  which  to  find  risetime') 

pause ( 2 ) 

% 

%  Determine  peak  absolute  value 

np  =  ne-ns+1; 
peak  =  0.0; 
for  I=ns:ne 
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if  abs(s(I,2))  >  abs(peak) 
peak=s (1,2) ; 
npk=I ; 
tp=s (I, 1) ; 

end 

end 

% 

%  COMPUTE  RISETIME  OF  PEAK 
% 

ys  =  abs(peak)  *  0.1; 
ye  =  abs(peak)  *  0.9; 

nye=]nax  (f  ind(abs  (s  (ns  :ne,  2  )  )  <=  ye  &  s(ns:ne,l)  <  tp))  +  (ns-l); 

nys=max (f ind(abs (s (ns :ne, 2 ) )  <=  ys  &  s(ns:ne,l)  <  tp))+(ns-l); 

dy=abs ( s (nye, 2 ) -s (nys , 2 ) ) ; 

dx=s (nye, 1) -s (nys, 1) ; 

dyr=ye-ys ; 

dxr =dyr /dy*dx; 

npr=nye-nys ; 

% 

plot (s(;,l),s(:,2),s (nys : nye, 1) , s (nys : nye, 2 ) , tp,peak, ' o' ) 

title ([ 'Risetime  is  ' ,num2str (dxr) , '  with  ' ,nuin2str (npr) , 'points '] ) 
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